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Abstract 

Aggregation of ice crystals is a key process gov- 
erning precipitation. Individual ice crystals exhibit 
considerable diversity of shape, and a wide range 
of physical processes could influence their aggrega- 
tion; despite this we show that a simple computer 
model captures key features of aggregate shape 
and size distribution reported recently from Cirrus 
clouds. The results prompt a new way to plot the 
experimental size distributions leading to remark- 
ably good dynamical scaling. That scaling inde- 
pendently confirms that there is a single dominant 
aggregation mechanism at play albeit our model 
(based on undeflected trajectories to contact) does 
not capture its form exactly. 



1 Introduction 

It has recently become possible to collect large sam- 
ples of high resolution cloud particle images in real 
time, opening up the modelling of cloud dynamics 
to detailed comparison with nature. Figure 1 shows 
ice crystal aggregates from a Cirrus cloud over the 
USA, captured by non-contact aircraft-based imag- 
ing. Such aggregates can be seen to be comprised of 
varied rosette ice crystal types, and detailed statis- 
tics have recently been published on both the clus- 
ter aspect ratios [Korolev and Isaac 2003] and size 
distributions [Field and Heymsfield 2003] in Cir- 
rus clouds. Such aggregation is a key feature of 
cloud development in the troposphere and can be 
quite crucial to the development of precipitation, 
whether it reaches the ground as snow or melts first 
to arrive as rain. 

The openness of the aggregates significantly ac- 
celerates their growth. Two clusters (labelled by 
pass approach with their centres closer than 
the sum of their radii r» + r, at a rate proportional 



to 



7r(n + rj) \vi 



(1) 



where for each cluster the sedimentation speed v 
varies inversely with its radius r and mass m as 



V/Pp f mgp \ 
r \ rj 2 J 



(2) 



Here r\ and p are the viscosity and density of the 
air, g is the acceleration due to gravity, and we 
have assumed that only one geometrical radius is 
relevant: Mitchell [1996] discusses an elaboration. 
Given the above, the rates of aggregation per unit 
time for fixed cluster masses vary linearly over- 
all with the cluster radii, and openness of aggre- 
gate structure enhances aggregation rates despite 
lowering fall speed. For real aggregates this is a 
significant factor: using data from Heymsfield et 
at, [2002], one finds that rosette aggregates 2mm 
across (which yield 0.5 mm droplets) aggregate four 
times faster than when melted. For cloud particles 
it is also relevant to consider the rates of aggre- 
gation per unit of distance fallen (rather than per 
unit time), which at fixed mass is proportional to 
the square of radius, leading to 16 times enhance- 
ment for ice over water in the example cited. 

We have made computer simulations of ice ag- 
gregation based on equations and (J2J, tracing 
trajectories through possible collisions to obtain ac- 
curate collision geometries. We assumed that all 
collisions led to rigid irreversible joining of clus- 
ters, as the openness of the experimentally observed 
clusters suggests little large scale consolidation of 
structure upon aggregation, and that cluster orien- 
tations were randomised in between collisions but 
did not change significantly during them. We took 
the sedimentation speeds to be governed by iner- 
tial flow, for which the mass dependence function 
in equation J2J is given byF(A) cx X 1 ' 2 . Details of 
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implementation are given in a longer paper [ West- 
brook et ai, 2004]. 

Some representative computer aggregates are 
shown in figure 1 alongside the experimental 
ones. Our simulations used three dimensional cross 
shapes for the initial particles as a crude represen- 
tation of the experimental bullet rosettes. 

Figure 2 shows a quantitative comparison of ag- 
gregate geometry, in terms of the ratio of cluster 
spans perpendicular to and along the direction of 
maximal span, as measured from projections in an 
arbitrary plane of observation. We find that differ- 
ent initial particle geometries (rosettes, rods) ap- 
proach a common asymptotic average cluster value. 
The aspect ratio of CPI images have been similarly 
calculated [Korolev and Isaac 2003], the results of 
which have been overlayed onto figure 2, and these 
appear to approach the same value. This universal- 
ity of aspect ratios provides direct support for our 
hypothesis of rigid cluster joining upon contact. 

A deeper indicator of universality is provided by 
the fractal scaling of ice crystal aggregates, where 
one tests the relation m oc £ df between aggregate 
mass m and linear span £. Our simulations and ex- 
perimental observations [Heymsfield et at, 2002], 
rather accurately agree on the fractal dimension 
d f = 2.05 ± 0.1 and in Westbrook et al. [2004] we 
discuss theoretical arguments leading to df = 2. 

Our simulations conform well to dynamical scal- 
ing of the cluster size distribution. This means that 
number of clusters per unit mass varies with mass 
and time together as a single scaling function, 

where S(t) is the weight average cluster mass. This 
relationship is confirmed in figure 3, where we 
rescale the mass distribution from different times 
in the simulations onto a universal curve. 

The scaling function which we observe in figure 
3 exhibits power law behaviour with &(X) oc X~ T 
for I < 1 with t w 1.6. This is not intrinsi- 
cally surprising (and indeed it matches theoretical 
expectations [Van Dongen and Ernst, 1985]) but 
it has forced us to abandon the way experimen- 
tally observed distributions of cluster linear size 
have hitherto been plotted. The problem is that 
given equation J2J) and its observed power law form, 
we must expect that the distribution of clusters 



by linear span t should at small £ take the form 
^(£,t) cx which diverges as £~ 2 - 2 us- 

ing our observed exponents. For small enough crys- 
tal sizes this behaviour will be modified by the role 
of growth processes other than aggregation, but 
that lies outside the scaling regime. 

Because of the divergence one has to take great 
care in constructing a characteristic linear size L(t), 

i k+i 

where the natural choices are L k (i) — 4^'" ater " 

and the lowest whole number k for which the de- 
nominator is not dominated by the smallest clusters 
is k — 2. The simplest natural scaling ansatz for 
the cluster span distribution is then found [ West- 
brook et ai, 2004] to be 

where M k (t) = J2dusters ^ ■ Figure 4 shows that 
this scaling ansatz works acceptably for our simu- 
lation data and well for the experimental observa- 
tions. The latter are rich data because cluster span 
is one of the simplest automated measurements to 
take. 

The rescaled distributions from simulation and 
experiment agree fairly well but not perfectly, as 
shown in Figure 4. One experimental reservation is 
the fall-off of experimental observation efficiency at 
small sizes, where clusters can be missed. However 
our scaling procedure itself is in effect expressly de- 
signed to avoid sensitivity to this, and the super- 
position of the experimental data down to small 
reduced sizes looks good. Indeed it looks so good 
that the transient flattening around £/L w 1, which 
is absent from the simulations, appears to be sig- 
nificant. 

One suggestion for the flattening around the mid- 
dle of the rescaled distribution is that it might be 
associated with the peculiar feature of the collision 
rate being zero between clusters of equal sedimen- 
tation speed. Our simulations include this feature 
but for low relative approach speeds, where each 
cluster has more time to adjust its momentum in 
response to the other perturbing the local airflow, 
our approximation of ignoring hydrodynamic inter- 
actions (and hence phenomena such as wake cap- 
ture) is less accurate. 

In summary, we have a fairly complete under- 
standing of the geometry of the atmospheric ice 
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crystal aggregates, dominated by sticking upon en- 
counter. Further details of the sticking mechanism 
(which we did not include) appear not to be impor- 
tant for the cluster geometry and the excellent scal- 
ing superposition of the experimental cluster size 
distributions suggests sticking efficiency does not 
favour aggregation at particular sizes. The sim- 
plest interpretation of these observations is that al- 
though the sticking probability might be low for a 
single microscopic event, many such contacts will 
be attempted during a cluster-cluster encounter so 
that eventual escape is unlikely. The actual sticking 
mechanism between ice crystals remains an intrigu- 
ing open question, particularly for the low temper- 
atures of figure 1. 

The fact that the same evolution is seen for differ- 
ing initial monomer populations (rods and rosettes) 
suggests that a single set of geometric relationships 
for ice aggregates can successfully be applied in a 
wide range of cloud conditions. This would lead 
to greater accuracy in retrieving cloud properties 
such as precipitation rate and predicting the radia- 
tive affect of ice crystal aggregates upon the climate 
system. 






* 




Figure 1: (a) Ice crystal aggregates images ob- 
tained from an aircraft flight through cirrus cloud, 
at temperatures from — 44°C to — 47°C (~9 km al- 
titude), using a cloud particle imager (CPI, SPEC 
Inc., USA). The pictures shown are aggregates of 
rosette ice crystal types, (b) Aggregates as simu- 
lated by our computer model which assumed rigid 
joining when clusters collide under differential sed- 
imentation. 
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Figure 2: Mean aspect ratio for projected ice 
agggregate images, where the aspect ratio is mea- 
sured as the longest span L max divided into the 
span perpendicular to the longest L perp . Grey lines 
show cloud data of Korolev and Isaac [2003] plot- 
ted against longest span in microns for a range of 
temperatures between 0°C and — 40°C. Black lines 
show simulation data plotted against longest span 
in arbitrarily scaled units, where the initial parti- 
cles were three dimensional crosses (solid line) and 
simple rods (dashed). 




Figure 3: Dynamical scaling of the cluster mass 
distribution from simulations. The superposition 
of data from different times supports equation J3J|, 
and the linear portion indicates $(X) oc X~ T at 
small X with exponent r = 1.6. 
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Figure 4: Cluster length distribution, rescaled as 
discussed in the text. The gray lines show experi- 
mental distributions at altitudes of 9.5km (— 50°C) 
to 6.6km (— 28°C) in the cirrus cloud of Field and 
Heymsfield [2003] obtained during an ARM (At- 
mospheric Radiation Measurement program) flight 
(9th March 2000). Each experimental size distri- 
bution represents an in-cloud average over 15 km. 
Black lines show simulation data. 



5 



